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We reconsider the treatment of Lise and Jensen (Phys. Rev. Lett. 76, 2326 (1996)) on the 
random neighbor Olami-Feder-Christensen stik-slip model, and examine the strong dependence of 
the results on the approximations used for the distribution of states p{E). 
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The work of Olami, Feder and Christensen |l[] on a slip-stick earthquake model indicated, some time ago, that self- 
organized criticality may occur without a local conservation law. Recently, it has been claimed by Lise and Jensen 
that the random-neighbor version of the OFC model also presents critical behavior above some critical dissipation 
level ac < ctQ, where ao — l/q is associated with local conservation for q neighbors. These authors based their claims 
on some theoretical mean field arguments and on numerical simulations with systems with up to A'^ = 400^ sites. In 
order to perform the mean-field calculation they had to make many different assumptions about the behavior of the 
model. 

More recently, Chabanol and Hakin |^ , and Broker and Grassberger |4| performed a more detailed analysis of the 
same model, showing that what has been interpreted as a critical behavior in |^ indeed corresponds to a subcritical 
region with very large (but finite) mean avalanche sizes. Although Broker and Grassberger Q| gives a comprehensive 
treatment of the random-neighbor version of the OFC model (which we will designate R-OFC), it may be of interest 
to detect exactly where the theoretical arguments given in |^ fail, since that point is not transparent in their paper 
and similar problems may occur or be of interest in the future. This is the aim of our paper. We will show that the 
problem is not in the method used in |^ (which eventually can give useful informations about the mechanism behind 
SOC) but in the strong dependence of the output of the calculations on the exact form of the distribution of states 
p{E) of the system. 

To reinforce the strong dependence of the results on the specific form of p{E), we revisit the R-OFC model, but 
this time introducing a simple and small modification on the p{E) distribution, that consists in replacing the interval 
[0, Ec] where the uniform distribution used by Lise and Jensen was defined by the interval [0, E*], with E* < Ec (Ec 
is the threshold value above which the sites become unstable and relax), that is 

Pi^) = ^ - (1) 

where Q{x) is the Heaviside function (see Figure (1-a)). 

The random version of the OFC model (R-OFC) consists of N sites initially with an energy Ei < Ec, for i = 1, N. 
The sites with energy E bellow Ec are stable sites (inative) and will be labelled by the superscript — ; the sites with 
energy E above Ec are unstable (active) and will be labelled by the superscript The energies of all sites are 
increased slowly until the instant t when the energy of a certain site i reaches the value Ec- This site becomes then 
unstable and the system relaxes in a very short time scale according to the following rules: 

E,{t + 1)^0, 

Emit + I) Emit) + aE+ , (2) 

where Em stands for the energy of q other sites chosen at random, and a < l/q. Eventually some of these q sites 
may become unstable, relax also, and so may generate an avalanching process that only stops when the energies of 
all sites are again bellow Ec- If we have a = ac = l/q we say that the system is conservative. 

Following B] , the probability of an inactive site to be activated by receiving a contribution aE^ of an active site is 
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The branching ratio a is the average number of new unstable sites created by a unstable site that relaxes. Clearly, 
in order to have a critical branching process, we must have (7 = 1. For q random-neighbors, we have 



a = q 
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where p{E^) is the distribution of states for unstable sites, that is, sites with energy above E* . Adopting the notation 
(. . .) = c , for averages, we fand 

^ ' J p(E) dE ' 
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The critical branching ratio a — 1 defines a value etc above which infinite avalanches may occur. The quantity {E^ 
is estimated H as 



1 — a 



The average on E is calculated as 
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E* -E^ + a {E+) 
With this result and Eq. (^), we get 

(2a - a^)x^ + [2£;*(1 -a)- 2Ec] x + E^ - {E*f = 0, 
where x — {E^). This leads to the simple solution 

Ec + E* 



{E+)^ 



2- a 



Finally, we get for the branching ratio 



a = q 



E*-E, ^ a{E, + E*) 
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The critical condition (7 = 1 leads to 
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Ec + E*l{2q) 

For example, li q — A (case studied in j^), we have 

_Ea- iE^/A 
E, + E*I8 ■ 
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(12) 



If we consider E* = E^ we recover the Lise and Jensen value ac — 2/9. But the value of ac has a strong dependence 
on the value of E* . We see that, already for the value E* /E^ = 24/25 = 0.96, etc achieves the physical limit 0.25. 

In Figure 2 we show the behavior of the normalized coupling die — qotcjA as a function of E* . With this normal- 
ization, the conservative case always corresponds to a = 0.25 (as for (7 = 4). We see that the allowed region of values 
for E* so that cic < 0.25 is very narrow for any g, and that the value of etc varies strongly in this region. We have 
already shown in another work (see |^) that already with lattices with N = 600^ sites it is possible to see a finite 
mean size avalanche for a = 0.23, contrasting with the results of based in simulations in lattices with N = 400^. 

Besides showing that Lise and Jensen's approach is not robust with respect to p{E), we may ask about what kind 
of model produces the uniform distribution 

p{E) ^ eiE)Q {1 - E) (13) 

used by those authors. We found that an extremal version of the Feder and Feder model, hereforth called EFF model, 
indeed produces this distribution. Extremal here is used in the same sense it was first used in the Bak-Snappen 
co-evolution model ||^. In the dynamics of extremal models there is no driving step. We locate the site i with the 
largest value of Ei — max{£'j} at the initial instant t. This site relaxes following the original FF rules: 

Er{t+l)^V, (14) 

Enn{t+l) =Enn{t)+a, 

where 77 is a noise, 1] € [0, e]. If we consider e = Q: = l/4we will have (p^. In this model the size of an avalanche 
is defined as the number of sites with energy iJ™ > 1 that relaxes in a row. Now, if we repeat the Lise and Jensen 
calculation using these EFF rules instead of the OFC rules, we get the self-consistent result ac = E^jq. 

It is also possible to show that a more realistic assumption about p{E^ leads to essentially the same results obtained 
by H,^. If we simulate the R-OFC model with q = 4 we will get a energy distribution p(i?) with four peaks They 
show clearly that "piE) is not a simple constant. We then decided to repeat the same calculations but supposing this 
time that p(Ef) had the (more realistic) form shown in figure (1-b), where Ap is half the width of each peak and Af, 
is the width of the gaps between two peaks. That means 

. f a, for E Ix OY E ^ I2 OY E ^ l-i or E ^ I4 , , 

^ 1 0, otherwise ^ ' 

where h = [0, Ap], h = [Ap+ A,,, 3Ap + Af,], I3 = [3Ap +2 A,,, 5Ap + 2A,,], and h = [5Ap + 3 At, 7\ + 3 A,,]. 

We also have that E* = 3Ab + 7Ap is the maximum value for which p{E) 7^ 0. 
Then we have 

/ P{E)dE 

P+ (E^) - ^ (16) 
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The lower limit of the integral in the numerator, Ec — aE^ , can now belong to any of the four intervals that define 
the peaks of the distribution, to which we will asssign the indices i = 1, 2, 3, 4. Considering each one of the possibilities, 
the integrals P]_ (i?^) assume the generic form 



The branching rate is given by 



cr = 4Fi = 4 
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In a similar way used to obtain eq. (|7|), we calculate an expression for {E ), which is associated with eq. (||) to 
lead 

[7Ap + (^-l)Aj,] (l~a) y/y- 
\ / a(2-a) a {2 -a) 2a (2 - a) ' ^ ' 



where 



y^=A {Ec (1 - a) - [7Ap + (i - I) At] f + 4a (2 - a) [x, - 14 (z - 1)] Ap A^, (20) 



with Xi = 24, 26, 32 and 42, for i = 1, 2, 3 and 4, respectively. Imposing the branching condition cr = 1 and using eq. 
(|l|) we get 

7Ap (2 + a) + 4 (i - 1) Af, - AE,, (1 - a) + TaE^ ± /l^ = (21) 

For instance, if we take Ap — 0.08 and A;, — 0.1 (that corresponds to Figure (1-b)), the critical branching condition 
leads to values of a* outside the physical range (that is, a* > 1/4). Therefore, in this particular example, it is 
physically forbidden to assume that cr = 1, so there is no self-organized critical state. 

If we take the limit for the conservative case (that is, Ap and Ab ctEc), the four peaks tend to four delta 
functions at {i — 1) a Ec, and it is easy to see that the condition a = 1 leads to the only possibility a* = ac = 1/4 
(we obtain a* > 1/4 for i — 1,2, 3). It can also be shown that, if we consider the limit Af, ^ and Ap — > Ec / 7, 
(that is p{E) is constant in the interval [0, Ec] which corresponds to the approximation of ref. |Q) then a* = 2/9. 

In general, for all values of i, the regions of the parameter space associated with a < 1/4 are determined by 

Ec - — A„ - =^Ab < 0. (22) 
24 ^ 21 ' 

From this inequality plus the relation Ec > 7Ap + 3 A;,, (see Figure 3) we see that only for a very small range of the 
parameters Ap and A^ there are values of ac in the physical range (0 < ac < 1/4). In all of those cases. A;, is very 
small and the shape of p{E) is very close to the constant form used by Lise and Jensen. Moreover, a varies strongly 
in these allowed ranges. 

In conclusion, we showed that, besides not having considered lattices big enough, the problem with the approach 
used by Lise and Jensen in is not in the method itself, but in the strong dependence of the output of the calculations 
on the compatibility between the distribution p(E) and the assumed dynamical rules which presumably lead to it. 
We also showed that the p{E) approximation used by Lise and Jensen does not correspond to the model they intend 
to analyse, namely R-OFC, but to another non-conservative model, designated by us EFF (Extremal Feder and Feder 
model). If we adopt the EFF's dynamical rules, the Lise and Jensen's method will lead to the right conclusions. In 
the end, we followed the same approach but now considering a more realistic approach for »(£^) and get essentially 
the same results that had already been obtained through the use of other arguments in ||,|[, that is in the R-OFC 
model there is SOC only in the conservative limit. 
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Figure Captions 

Fig. 1- a) The uniform approximation for the energy distribution p{E) of the R-OFC model. The soUd hne 
corresponds to the interval of E used by Lise and Jensen, [0, Ec]; the dashed line corresponds to a smaller interval of 
E, [0,E*], used in our calculations, b) A more realistic approximation (non-uniform) of the energy distribution of 
the R-OFC model, for g = 4 (with four peaks). The parameters Ap and are, respectively, the half- width of the 
peaks and the width of the gaps. The value of p{E) at the peaks is a and Ec > E* = 7Ap + 3Ai,. 

Fig. 2- Normalized 'critical' dissipation level etc = ka/A as a function of E* . Values for ac above 0.25 are not 
physically admissible. 

Fig. 3- Space of parameters for p{E) in terms of 7p — Ap/ Ec and 7b = A-hjEc- The shaded regions correspond to 
the intersection between a < 1/4 and (lAp -|- SA^) jEc = l^p -|- 2>^h < 1- Depending on the value of Ec — aE~^ , we 
have (a) Ec - aE+ e [0, Ap]; (b) Ec - aE+ e [Ap + A^, 3Ap + A^]; (c) Ec - aE+ e [3Ap -|- 2Ab, 5Ap + 2Ab]; and (d) 
Ec - aE+ e [5Ap + 3Ab, 7Ap + 3Ab]. - 
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